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Abstract: 

The appearence of long-range correlations near the Dirac point of a Dirac-like spinor model with 
random vector potential is studied. These correlations originate from a spontaneously broken symmetry 
and their corresponding Goldstone modes. Using a strong-disorder expansion, correlation functions and 
matrix elements arc analyzed and compared with results from a weak-disorder expansion. The local 
density of states correlation and the overlap between states above and below the Dirac point are charac- 
terized by a long-range behavior. The correlation range decreases with the distance from the Dirac point. 
Transport is diffusive and the diffusion coefficient is proportional to the one-particle scattering time for 
any strength of disorder. A consequence of the special properties of particle-hole scattering is a constant 
microwave conductivity for weak as well as for strong disorder, describing a deviation from conventional 
Drude-like transport. Some properties of the model can be linked to a kind of Kondo scale, which is 
generated by disorder. Finally, the properties of the wave functions at the Dirac point are characterized 
by their participation ratios, indicating a critical state at the Dirac point. 

PACS numbers: 81.05.Uw, 71.55.Ak, 72.10.Bg, 73.20.Jc 

1 Introduction 

The electronic properties of graphene are closely related to the existence of a spinor wave function and 
the two Dirac nodes in the band structure [U [2J [3] . This implies an unconventional behavior which is 
associated with the Klein paradox [4]. An important effect is the scattering between the hole and the 
particle sector of the Dirac cones, leading, for instance, to the zitterbewegung El [10] . A constant 
contribution to the microwave conductivity is another consequence of this effect, causing a deviation from 
conventional Drude-like transport [TTJ [TH [T5] . The latter has been studied only for weak disorder. It 
would be interesting to analyze these effects, their modification or even their destruction in the presence 
of strong disorder. 

The presence of disorder in graphene has been discussed in a number of recent papers, considering 
different physical conditions. It seems that disorder should appear in the effective Dirac-like Hamiltonian 
in the form of a random vector potential due to instability of the translational order in 2D [Til 151 [T51 [TBI [17] . 
Scattering by a random vector potential affects only the phase of the wave functions. Its influence on 
transport is such that no localization has been observed, even for very strong disorder [18j . However, it is 
believed that disorder has a substantial effect on the magnetoresistance, observed as a suppression of the 
weak-localization peak [T4J, [19] . Another effect, which might be caused by disorder, are fluctuations of the 
charge distribution near the Dirac point. Recent experiments have revealed that in graphene long-range 
charge fluctuations exist [2D] . 

The Dirac point is very special because it is associated with zero-energy modes and a vanishing density 
of states. Consequently, it might be very sensitive to disorder. It is known that, depending on the type 
of disorder, the properties of Dirac fermions are strongly affected by randomness |21| I22j. An example 
is the prediction of a dramatic change of the density of states due to disorder. On the other hand, 
some quantities, like the conductivity, are remarkably robust [531 1211 HH US] . Of particular interest is 
the question how disorder affects the interband scattering. Weak disorder has no effect on particle-hole 
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scattering, only on the inter-particle scattering (cf. [H]). All this indicates that the Dirac point is an 
interesting object to study in more detail, especially when disorder is strong. 

To understand the effect of a random vector potential in a Dirac-like Hamiltonian, it has to be noticed 
that the Dirac point is characterized by a chiral symmetry. We discuss this symmetry and its spontaneous 
breaking by the random vector potential. It is known that the spontaneously broken symmetry implies a 
robust diffusion mode for weak disorder [24] and for strong disorder [18j . In this paper we will demonstrate 
that the local density of states correlations have a long-ranged behavior, although the average density of 
states is finite. Some properties can be linked to a kind of Kondo scale, generated by disorder [26l[27]. The 
analysis of the scattering between different parts of the Dirac cone leads to correlations between different 
energy states, whose correlation are expressed in terms of transition matrix elements. A special case 
is the corresponding matrix elements of the position operator with respect to ±E energy states, which 
appear as important contributions to the microwave conductivity. Finally, the second moment of the 
local density of states can be associated with the inverse participation ratio. This quantity characterizes 
the properties of the wave functions at the Dirac point. 

In order to study the effect of symmetry breaking in the presence of disorder, we apply a perturbation 
theory for strong disorder. Strong-disorder expansions have been quite successful in different fields of 
physics. In the case of Dirac fermions it has been observed that the expansion has a number of advantages 
in comparison with a weak-disorder expansion [18] . The main reason is that the leading order of the 
expansion is quite simple. Nevertheless, the strong-disorder expansion cannot be applied directly to the 
one-particle Green's function of the Dirac-like Hamiltonian because it has, in contrast to a Schrodinger- 
like Hamiltonian, divergent contributions from self-crossing loops. (The strong-disorder expansion can be 
formulated as an expansion in terms of the hopping elements of the Hamiltonian.) These problems can 
be avoided by applying the expansion to the saddle-point integration. Then the integration is restricted 
to the most important contributions, namely the saddle-point manifold, that takes into account the 
underlying symmetries of the model. In this case the strong-coupling expansion can be controlled and 
leads in leading order to a model with diffusion propagator. It describes the Goldstone modes which are 
generated by spontaneous symmetry breaking. The results of this expansion can be compared with those 
of the more common weak-disorder expansion. The latter gives in leading order a factorization of the 
average two-particle Green's function into a product of two average one-particle Green's functions, which 
can be treated within the self-consistent Born approximation. Both approaches are compared to study 
similarities and differences of both regimes. 

The paper is organized as follows. In Sect. 2 we introduce the generalized Dirac model and discuss 
spontaneous symmetry breaking. This includes the definition of some physical quantities, which are 
studied in the rest of the paper. The framework of functional integration for our model is presented in 
Sect. 3, and the corresponding symmetries are analyzed in Sec. 4. Then the averaging with respect to 
the random vector potential and the related saddle-point integration are performed (Sect. 5). The results 
of the calculations are discussed and compared with weak-disorder approximations in Sect. 6. 



2 Model and spontaneous symmetry breaking 

In this paper we study the generalized Dirac Hamiltonian for a spinor-1/2 state in a random vector 
potential 

H = Ho +vai, H = hicri + h 2 <T 2 ■ (1) 

a j (j — Oj 1)2,3) are Pauli matrices, hj is an operator that acts in space, and v is space-diagonal with 
random variables v r . The latter are independently Gaussian distributed with zero mean and variance g. 
Special cases for Hq are 2D Dirac fermions, whose Fourier components are 

hj = kj 

and the tight-binding model on a honeycomb lattice with 

3 3 

hi = — t y cos(aj • k), hi = —t/ sin(aj • k) 



with the lattice vectors of the honeycomb lattice 

ffi = (->/3/2,l/2), a 2 = (0,-l), a 3 = (\/3/2, 1/2). 

i/o can be diagonalized as i?o = diag(ek, — e^) with efc = \J h\ + h\. The Hamiltonian iJ is invariant 
under the continuous transformation _ff — * e la<73 H e IQCT3 . This symmetry can be broken spontaneously. 
To measure spontaneous symmetry breaking (SSB), we consider the one-particle Green's function 

G(z) = (H + z)- 1 

and calculate the difference of this expression before and after the symmetry transformation. For the 
special case with parameter value a = iir/2 this reads 

GO) - ia 3 G{z)ia 3 = (H + z)" 1 - (H - z)' 1 . 

Since z breaks the symmetry, we send it to zero. If we choose z = ie, the results for the diagonal elements 
of G are proportional to the local density of states at the Dirac point: 

p r (0) = S rr (H) oc \im[{H + ie)" 1 - (H- ie)- 1 ] . 

Thus, a nonzero p r (0) indicates SSB. This is similar to spontaneous symmetry breaking in classical 
systems, e.g., in the case of the magnetic phase transition of a classical Heisenberg model [28j. However, 
there is a remarkable difference: In the Heisenberg model there is SSB only for systems with dimensionality 
d higher than 2, whereas this restriction does not hold in a disordered system. Even for d = 0, where 
H = vcri, we have SSB, since the DOS reads 

(p(0)) = lim / 6 P(v)dv = lim / - 1 ^ PUx)dx > 
^ ;/ e-,0 J e 2 +v 2 w moJ 1+i 2 v 1 

for a continuous distribution function P(v), e.g. a Gaussian. Another interesting point is that the second 
moment of the density of states diverges like e~ 1 : 

(P(°) 2 ) = / {e 2+ v 2 ? P (v)dv = e- 1 J -^-L^P(ex)dx . (2) 
This simply reflects the fact that the distribution of the density of states has a Lorenzian form. 



2.1 Physical quantities 

A characteristic property of the Hamiltonian H in Eq. ([1]) is that its spectrum is symmetric with respect 
to the energy E = (the Dirac point), consisting of states ^±e in the upper and the lower band (particle 
and hole parts of the Dirac cones). These states are connected by the transformation ^-e — 03 '5 e 
because the eigenvalue equation 

H^e = E^ E 

gives 

a3Ha3a 3 ^E = Eaz^ e 

which becomes with a^Ha 3 = — H 

Ha 3 ^ E = -Ea 3 y E ■ 

Scattering between ^e and ^-e leads to interesting physical phenomena, like the zitterbewegung [51(5] or 
a constant contribution to the microwave conductivity [111 1291 1 1 3] . The latter can be seen by considering 
a typical contribution to the frequency-dependent conductivity at frequency uj (cf. Appendix A) as 

2 rU>/2 

= -7T W / ? r 2 k Tr 2 [G r0 (uj/2 -E- ie)G r(-u/2 - E + ie) 

4 « J-u/2 r 

+ G r0 (-uj/2 + E- ie)G 0r {uj/2 + E + ie)]dE , (3) 
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where Tr 2 is the trace with respect to Pauli matrices. This can also be expressed in terms of the matrix 
elements of r\ between the states <&± u /2-e (cf. Appendix B) as 



e 2 r /2 

(jII ^^~Jj; UJ (( $ --/2+Bkfcl^/2+£> + (^/2-£kfcl$-c/2-£»^ • (4) 

4/1 J-u/2 



For a system at the Dirac point E = 0, we expect that the main contribution to the conductivity comes 
from the scattering between $_ w /2 and < & w /2, where the scattering does not change momentum (cf. Fig. 
1). Therefore, we consider in the following only 

2 2 

*o(w) = -|^w 2 <$_ L , / 2|^.|$ L , /2 ) + ($ w / 2 |rg|$- w /2» = -^ 2 ($- u /2\rl\$ u / 2 ) , (5) 

where the last equation is a consequence of the symmetry of the matrix element. Then a constant 
conductivity, found theoretically 112] and experimentally [13] , requires a matrix element that diverges 
for oj ~ as 

(*-«/a|r2|«a;/2) ~^~ 2 . 

There are two other correlation functions that are of interest for the characterization of the disordered 
system of quasiparticles in graphene. One is 

Crr' = (Tr 2 [5 rr ,(H ~ E)5 r , r {H - E))) (6) 

which is related to the DC conductivity at T = to — through the expression [25] 

*w = 2tt^ lime 2 £(r M - r'^C rr> . (7) 

r' 

The other correlation function is 

D rr , = (Tr 2 [6 rr (H - E)]Tr 2 [8 r , r ,(H - E)}) (8) 

which is the correlation of the local density of states at sites r and r'. For r' = r, D rr is related to the 
participation ratio and the inverse participation ratio (cf. Appendix C). These are quantities that measure 
the statistical properties of the wavefunctions and their localization behavior [30j . The participation ratio 

p (2) _ (Pr) 2 



(P 2 r) 



vanishes for localized states because the local density of states p r has a broad distribution with divergent 
second moments in this case, similar to the single-site density of states of Eq. |(5J). Using D rr the 
participation ratio reads 

P < 2 ) = lim M . (9) 

e— »0 L) rr 

Moreover, the inverse participation ratio is also related to the second moment of the local density of 
states (cf. Appendix C): 

P (2) = lim eD rr . (10) 

It vanishes as one approaches the regime of extended states, coming from the localized regime. 
The density operator 5(H) can be expressed by the one-particle Green's function 



G ± = (H±ie)- 1 

as ^ 

S rr '(H) = — (G rr , - G rr ,) . 

With the relation a^G^cr^ = —G T the delta function also reads 

4 „., (ff) = i [GX ,+. sC>sl = i(^.n » 
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This relation implies for the correlation functions at the Dirac point E = 

Crr ' = ~~2 ( G ?r',jj G t'r,jj) ' ( U ) 



7T^ 

4=1,2 

Orr' = -\{Tr 2 [G+]Tr 2 [G+ r ,]) . (12) 
Moreover, away from the Dirac point at energies E — ±uj/2 it gives 

A rr ,(u) = -(Tr 2 [a s G rr/ (u)/2 + ie)a s G r , r (u>/2 + ie)]) , (13) 
which can be used to write for the conductivity of Eq. ([5]) 

2 2 

o-oH = -^^ 2 5I^^o(^) = -|^^ 2 (^ /2 |r2|$_ w/2 ) . (14) 

r 

3 Functional Integral 

All three correlation functions in Eqs. HJ), (jX3J) and (TiU]) are written as products of two Green's functions. 
We can express these products, before averaging, as a Gaussian functional integral with two independent 
Gaussian fields, a boson (complex) field Xrk and a fermion (Grassmann) field ^f r k (k — 1,2) and their 
conjugate counterparts Xrk and ^f r k- 



G rr >jf(ie)G r/rt k'k(it) = J Xr'fXrj^rk^r'k' e^p(- S (E))T>^T>x , (15) 
where Sq(E) is a quadratic form of the four-component field <j) r = (x-ri, Xri, ^ / r2, ^2) 

So(E) = -i^2(f> r -{H + ie + E) r , r 4 r ' (e > 0) . (16) 

The use of the mixed field <f> r has the advantage that an extra normalization factor for the integral is 
avoided. The extended Hamiltonian H = diag(H, H) of So acts in the boson and in the fermion sector 
separately. Using Eq. (|15p . the correlation functions now read 

Crr' = ~^Y1 [ Xr'jXrj^rj^r'j(exp(-So{0)))V^V X , 
3 



and 



D rr'=-\Ys I XriXrj*r'j'*r'i'(exp(-5o(0)))X>*X>X 
71 3,3' 



In the case of C rr ' and A rr >(oj) we arrange the products of the fields such that pairs at the same site are 
neighboring factors: 

Crr ' = ~2 J2[ Xr'j^r'3^rjXr3{exp(-S (Q)))V^Vx , 
3 

A rr ,{u) Xr'j'*r'i'*rjXrj(exp(-5 (a;/2)))OTOx • 
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4 Supersymmetry 



The Hamiltonian H in So is invariant under the transformation 

TT tj t \t t tt / ' n /e^+«> 3 (l + #/2) e^+P'^^a 3 , 
U = I7 (g,p)Usl7o(g ,*>) = ( g^'+p)^^ e «W)«(i - ^/2) ' ( 1 



U s = exp j , U (q,p) - 



with 

IF "~ ^ cn- : . a I' ' m, '" r ' - V o 
For 7j = diag(aj, <jj) there is the symmetry transformation 

u 7j u' = uu-S, = l3 

with 

U = U {q,p)V s U (q,p) = e ,( P + P > 3 (1 _ ^/ 2 ) J ■ ( 18 ) 

This implies the invariance 

H -> UHU' = H . (19) 
To simplify the notation one can introduce the phases tpij with 



and write 



fn=q + q', <pi2 = q + p', v>2i = q'+p, ^22 =p + p' 



e 2itpntT 3 Q _|_ _|_ e 2i^i2CT 3 ^/ J ^ ) Jgj(Vl2+^22)o'3 _|_ g*(Vl 1 + V21 )&3 j ^ 



<T3 



I r c «(¥'ll+V2l)o'3 _|_ g«(Vl2 + V22)CT 3 j^ ;(T3 g2j¥3220- 3 M _ _ g2iV2lO-3^,^, I ' 1 



5 Averaged correlation functions 

Averaging Eq. I|15p over the Gaussian distribution of v r means replacing exp(— So) by (exp(— So)) on the 
right-hand side of the equation. The latter can be written again as an exponential function (exp(— So)) = 
exp(— Si), where the new function Si contains also quartic terms of the field <f>: 

Si = -i^2<pr ■ (H + ie + E) r y(f> r > + g^ifir ■ li4> r ) 2 ■ (21) 

r,r' r 

Then it is convenient to transform the integration variables (Hubbard-Stratonovich transformation) as 

( XrXr Xr^r \ n _ / Q V &r \ f ™s 

where Q r , P r are symmetric 2x2 matrices and r , Q r are 2x2 matrices whose elements are independent 
Grassmann variables. Now the correlation functions can be rewritten as correlation functions in the new 
field Q r : 

°rr' = I (*i e )f JjfaQJrjj exp(-S 2 )2?[Q] (23) 

A rr ,{u) =-^2 J Tr 2 {a 3 ai& r/ )Tr 2 (a 3 aie r )exp(-S 2 )V[Q} (24) 
^ [ Tr 2 (-iaiP r ,)Tr 2 (aiQ r )eM~S 2 )V[Q] (25) 

TT J 



and 



o 
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with 



S 2 = V -Trg(Q2) + ln[detg[H + ie + E - 2 7l Q]] . (26) 
,9 

r,r' 

Trg is the graded trace, Tr2 the trace with respect to Pauli matrices, and detg the graded determinant 
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5.1 Saddle-point manifold 

The integration in Eqs. (f2"3")l , and (j2"5)) can be performed in saddle-point approximation. The saddle 
point (SP) is obtained as the solution of 6S2 = 0. Assuming a solution of the form 



Qo 



.V 
2 ' 



we obtain the parameter 77 from the SP equation 

7] = igTr2G rr (E + ie + iry) . 

For e = E = the SP equation is invariant under the global symmetry transformation 71 Qo 
of Eq. P^|) . This transformation leads to the SP manifold 

Qr = -i^lUrU^. = -*2 U 'r VlU, , 



(27) 

U 7l Q U' 

(28) 



where U r and UJ, are obtained from Eqs. (| 17[) and (|18p by replacing the transformation parameters if), 
... by space-dependent variables ip r , ... The form of QJ,, which is dictated by the symmetry of Sect. 4, 
implies for the action on the SP manifold that (i) the quadratic term in S2 vanishes and (ii) the remaining 
term becomes 



S' = lndetg(H + ie + E + irjUV) = lndetg(U -1 (Ho + ie + E)V'~ 



■irj) 



(29) 



This action contains the symmetry breaking field ie+E. Expansion in powers of r] 1 yields, after renaming 

u- x ->u, 

S' = ( —^Trg{\5\]') + ^Tr 9 (U'UH U'UH ) + 0{^) . 

Rescaling tp — > ij ll fi an d ^ — * V^ 1 '4 ) does not change the integration measure and allows us to perform 
an expansion of UU' in powers of r\~ x up to o(r]~ 3 ): 



UU' = 



(T + 2^ikt 3 /?7 - 2(pl 1 o- /r) 2 + 2ipipa /r) 2 2ipa 3 /ri + 2i(ip n + ^22)^0 /rf 

2ipa 3 /r] + 2i(ipn + ^22)^0 j 'if <t + 2^220-3/?? - 2(pl 2 a a /ri 2 - 2V>V>co/V 



(30) 

where we have used </5i2 + <^2i = q+p' + q'+p = <p>\\ + p>22- This provides an expansion of the action, 
where the leading order is a quadratic form in terms of the fields ip and tp: 

S' = ^ K rr'{t - iE)(-<pil tr ipil t r> + tP22,r<P22,r> + 2ty T $ r i) 



+ — K rr ,{Q)(ipu r + ¥?22,r)(</>ll,r' + f22,r')^rA' + Ofa 7 ) 



with 



K rr >(e - iE) = 4?7~ 



|^ ^ ^ hj mrr hj mrr I Z ^ Spy' ^ .?"?"' hij^T'y 

i=l,2 f j=l,2 



(31) 



(32) 



There is only one term that couples the Grassmann field ^ with the field y?. It turns out (cf. Appendix 
D) that this term drops out after the integration over tp. Moreover, the Jacobian J of the transformation 
should be J = 2~ N + o(r)~ 2 ) (N is the number of lattice sites) in order to satisfy the condition 



J e~ s ' JV[iP,ip] = 1 



(33) 



The symmetry-breaking term e — iE appears as a prefactor of a diagonal term. This action can now be 
used to calculate the correlation functions in Eqs. ([23]) . (|24l) . and ([25]) . An expansion of the components 
of the matrix field yields 

1 2 

o x Q = --[770-0 + 2itp n a 3 + -(-<£>n + ipi>)<ro] 

2 77 
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1 2 

-ia x P = --[wo + 2i<P22<?3 (V22 + V^Vo] 

2 r\ 

01Q = -i[a 3 + -{(fin + ^22)00)]^ 

(Ti6 = ~i[a 3 + -{(fin + ^22)00)]^ • 
V 

Integration over the fields tj) and if with respect to the quadratic action S' in 

(...)= J ,..e- s 'v[^ V ] (34) 

leads to (up to o(?7 -1 )) 

(criQ) = i<Joi]/2, -i(aiP) = icr a r]/2 , (35) 

since the integration rules imply 

+ (#) = (^ a ) + (#) = . 



It is also a consequence of Eq. (|33|) . Moreover, we need for the evaluation of the correlation functions in 
Eqs. ((23|) - ([25|) the following expressions: 

((a-l©r')jj(<7i<9r)^) = -(Vv^r) +o(rr X ) 
((^(o-ao-ie^JTra^ao-ier)) - -4<^ r ) + o^- 1 ) 

(Tr 2 (-i<TiPr>)Tr 2 {criQ r )) = -i] 2 - 4(VvVV') (VV'Vv) ■ 
2 , , T , 1 



and 

This leads to 
and 



C rr , = _^<^ r > = (36) 

A-r'M = -4r<VV</v) = 4*7r' ■ ( 37 ) 

9 9 



Finally, the correlation function of the local density of states reads 



D rr> = -^(j? + 4<VWV>(VV<M) = 4^2 + ^Kr-Kw ■ (38) 



In summary, all the correlations are expressed in terms of the inverse of the matrix K of Eq. 
Interesting is that only correlations of the Grassmann field ip appear, whereas the real fields ipjji do not 
contribute, at least in the approximation up to o(ry _1 ). 



6 Discussion 

All three quantities in Eqs. (f3T>l) . (I3T)) . and (f3"5)) are related to the same correlation function, namely to 
K~, — —2(ip r ?p r i) of Eq. (|32)l . The latter, or its Fourier components 1/K(q) with 

K(q) = -^[e' + c(q)}, e' = e-iE, (39) 

can be considered as the propagator of the average two-particle Green's function (G rr >G r ' r ). It describes 
the motion of two particles, created at the same time at site r' (cf. Appendix B). Although the two 
particles are independent (i.e., we ignore their Coulomb interaction here), the averaging over the random 
vector potential creates an effective interaction between them. This interaction is presented by the quartic 
term in S\ of Eq. (|21| . whose strength is g, the variance of the Gaussian distributed vector potential. A 
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consequence of the interaction is that the two-particle propagator describes diffusion, when we consider 
e' = iuj/2 and study q ~ 0: 

f 77 3 /4 



K(q) ioj + Dq 2 



with the diffusion coefficient 



D= l_d 2 c{q) 



2 da 



2 



q=a 2 V J k ^ [dh dh dkf 3 ) ■ [ ' 



Here it has been assumed that D is isotropic. 

Assuming weak disorder, the action S' in Eq. (j29|) can be expanded in powers of 77. (This weak- 
disorder approach is valid even for e = co = 0, in contrast to the factorization approach of Sect. 16. 5| A 
diffusion propagator was also found for this case, with a the different diffusion coefficient though [25 : 

A, = — 

for Dirac fermions and E = 0. Thus, the physics of the average two-particle Green's function is diffusive, 
both for weak (i.e. for 77 <§C 1) and for strong disorder (i.e. for r\ ^> f). The diffusion coefficient depends 
on the one-particle scattering rate r\ and the Hamiltonian Hq, given by its components hj in Eq. |T]). It 
is remarkable that both diffusion coefficients, the one of the strong-disorder expansion in Eq. ([40]) as well 
as D w of the weak-disorder expansion, are proportional to the one-particle scattering time r = i]^ 1 . The 
latter can be calculated from the self-consistent Born approximation [31]. However, D always decreases 
with increasing disorder, whereas D w is not monotonous with g but has a minimum due to the extra factor 
g. For realistic values of g, where the variance of the Gaussian distribution is < g < 1, the diffusion 
coefficient D w decreases with g. In any case, the increasing behavior of D w is beyond the validity of the 
weak-disorder approach. 

The integration with respect to the SP manifold of Sect. 5.1 allows us to identify the relevant 
transformation parameters for the long-range correlation functions. From the results in Eqs. (|36[) - (|3"8|) it 
is obvious that only ip, i.e. the transformation that mixes fermions and bosons, is relevant. On the other 
hand, the transformation inside the bosonic and inside the fermionic sector, provided by the parameters 
ipjji, is not relevant for the long-range correlations. This is in agreement with previous calculations for 
the conductivity, where the integration with respect to tpjji has not been taken into account [32] . If we 
project the symmetry transformation by choosing ip — t/j = m Eqs. (1I7|) and (fl8|) . there would also be 
a Goldstone mode, which becomes massless as we send the symmetry-breaking term to zero (i.e. e — > 0). 
This case has two interesting consequences: (i) the average density of states would be divergent at E = 
due to 

(0-iQ r ) = -iv<?o({<Pii, r ) - V 2 ) , 

and (ii) the long-range behavior of C rr i, A rr > and D rr i would disappear. 

For the calculation of the physical quantities of Sect. 5.1 we need the following expressions. The 
diagonal elements of K~ Y diverge logarithmically with e ~ as 

K - = Im ~ KoHeh 

which implies a divergency of the second moment of the local density of states D rr . Summation over the 
lattice sites gives 

I if 



K(q = 0) 4e' 

Finally, for the matrix element of r 2 , and the conductivity we need 

d 2 I rf 



(41) 



,=0 2e 
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6.1 One-particle scattering rate 77 

In the following the one-particle scattering rate 77 will be discussed for the specific case of Dirac fermions. 
Then the SP Eq. (J27J) reads 

A 2 



Zir 



1 



(43) 



where A is the momentum cutoff of the Dirac fermions. For weak disorder (i.e. g <C 1), the scattering 



rate is also weak. To study the Dirac point, we rewrite Eq. (J43f as 

n = -iE + . . (44) 

^exp(2nr]/g(iE + 77)) - 1 

and take the limit E = 

m = . 1=7= 1 , (45) 

This is shown in Fig. 2. For \E\ <C 770 we can solve the SP Eq. (|4~4"|) by an expansion in E: r\ = ?y + o(-E). 
On the other hand, for \r}\ we can iterate Eq. (|43[) with the initial value 

7/ = |-i£;in[l-A 2 /£ 2 ] = ^ ^ g n(£) + ^hi[A 2 /i? 2 -l]") (E 12 < A 2 ) . (46) 



77 is a reasonable approximation of 77 if 



l + -ln[A 2 /E 2 -l] 



7T 



< 1 



6.2 Density of states 

The density of states is proportional to the average one-particle Green's function 

p(E) = Tr 2 (G rr {E - it)) 

and describes two important features of our model, the spontaneous symmetry breaking and the one- 
particle scattering rate 77 in the SP approximation (cf. Eq. P7|) . In terms of the SP integration of the 
functional integral, Eq. (|35[) provides a finite average density of states at the Dirac point 

p(E - 0) - -i-T^atQ) = --Tr 2 {G X P) = VQ /g , 
9 9 

which is practically zero for a larger regime of g (cf. Fig. 2). Moreover, for \E\ ^S> \rj\ we get a linear 
behavior from Eq. (|46f 

p{E) « -Re{n) = J|i , 
9 2 

which reflects the density of states of pure Dirac fermions. Both results are in good agreement with a 
self-consistent calculation of the average density of states [27j [33] . 

Correlations of the local density of states D rr i have a long-range behavior. The corresponding Fourier 
transform 

fl(9)= l«W) (47) 

is a function of q 2 and diverges like q~ 2 at q = 0. It has has a kink (or a shoulder) at the edge of the 
one-particle spectrum A (cf. Fig. 3). 
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6.3 Microwave conductivity 

According to Eqs. HU), (|37[) and (|4^|) . the matrix element of r\ reads 

(3«/a|r2|3- w /a) = $>^ r0 H = -^4^ ■ 
Using the expression of the diffusion coefficient in Eq. ([4"0]) and taking e — * 0, the matrix element becomes 

<^i^/»>-^/.s:(^-9^)- ,48) 

Lower frequencies w (i.e. lower energies) are more important for scattering than states of higher energies 
due to their larger matrix elements. Moreover, if we interpret the matrix element as a measure of 
localization, the states with u> > are localized on a scale 1/ui. The matrix element, together with Eq. 
(fl~4|) . gives for the microwave conductivity 

»M " -^ 2 <*»/^l*-»/*> ~ t " |»') <« « »> ■ < 48 > 

This result indicates a constant microwave conductivity, at least for lu -C ?7o, since the prefactor cj 2 is 
compensated by the ui~ 2 behavior of the matrix element of r|. It should also be noticed that lim w _>o co(w) 
gives the DC conductivity of Eq. |7|) . This follows immediately from the definitions of the two conduc- 
tivities in Eqs. 0, dHJ) and from Eqs. ([36]), ([37|. 

The w -2 behavior of the r 2 matrix element in Eq. (|48[) obviously does not depend on the special form 
of hj, as long as the spinor structure of Hq exists. This implies that also for a parabolic kj dependence 
(e.g. in the case of a graphene bilayer [34], [35] [36] ) , the cancellation of the ui 2 terms in the conductivity 
takes place. 

6.4 Participation ratios 

The inverse participation ratio of Eq. pU]) vanishes like 

P (2 » ~ e(lne) 2 , (50) 
which indicates the existence of delocalized states at the Dirac point. The participation ratio 

P (2 Mlne)- 2 , (51) 

on the other hand, vanishes logarithmically. These two results are consistent with a critical point at 
E = 0, where there is a transition from localized to extended states. 

6.5 Perturbation theory for weak disorder and factorization 

The fact that diffusion is controlled by the one-particle scattering rate r\ raises the question about the 
quality of the one-particle approximation. The latter has been used frequently by factorizing the two- 
particle Green's function [37] [38] E] [15]. This approximation should be valid for weak disorder. It 
is based on the assumption that the Green's functions are uncorrelated and the averaged product is 
approximately the same as the product of the averaged one-particle Green's functions: 

This allows us to treat the average one-particle Green's functions within the self-consistent Born approx- 
imation: 

(G ± )^(H ±t V )- 1 = (H ±i V )- 1 , 

where r\ the imaginary part of the self-energy (or inverse scattering time) determined for Dirac fcrmions 
in Eq. (|4"3"]) . Consequently, the one-particle Green's function decays exponentially on the scale rf 1 . The 
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correlation function D rr / is constant and proportional to rj 2 , and there is no divergence for r' = r. This 
means that we have lost in the factorization the correlation term 



g 2 ir 2 



of Eq. ([55]) . These are the substantial differences between the strong-disorder expansion and weak- 
disorder perturbation theory at the Dirac point. If we go away from the Dirac point, we can study the 
matrix element rf,, approximated by the factorization as 

(^/ 3 |r2|4- w / 2 ) = - £ r 2 Tr 2 [a,G r0 {u/2 + ir,)a 3 G 0r (u;/2 + it,)] ~ { ~f\ ^ >> W , 
^— ' l_ — 4cj for ?7o "C w 



where we have assumed that a; <C A. In contrast to the result of Eq. (I48p . the matrix element does 
not diverge now if we approach the Dirac point w = 0, since there is the finite limit t/q 2 . This reflects 
the finite decay length 77^ 1 of the average one-particle Green's function. Away from the Dirac point the 
factorization works better. This can be seen if we insert the matrix element of r 2 into the conductivity 
ofEq. CHD 

2 2 2 

o^) = -^u 2 {<b^ l2 \rl\<S^ l2 )~^ = 2 e - (r, «^) (52) 

in the perturbative regime. Thus again, the conductivity does not depend on the frequency. Actually, 
the validity of the perturbative regime for the matrix element of r\ and the conductivity is big in terms 
of g due to 7j ~ e _7r / 9 . 



7 Conclusions 

Spinor states, described by a two-dimensional Dirac-like lattice Hamiltonian, were studied in an uncor- 
rected random vector potential. Our calculation, based on a strong-disorder expansion, has revealed 
that the quantum states develop long-range correlated fluctuations. In other words, the quantum sys- 
tem transforms the uncorrelated fluctuations of the random vector potential into long-range correlated 
fluctuations, for instance, of the density of states. The origin of this behavior is spontaneous symmetry 
breaking, which develops a massless (long-range) mode. The spontaneous symmetry breaking is measured 
by the one-particle scattering rate 77 or the density of states 77/ g. An important scale for this effect is 

which depends on the variance of the Gaussian fluctuations of the random vector potential g. (Although 
our approach is not valid for very small g (cf. Sect. 6.5), this scale is short for all reasonable values 
g rs 1.) r/Q separates regimes that are controlled by disorder (77 <C u>) from that which is controlled by 
energy (or frequency) (rj <C u>). For instance, a central quantity is the one-particle scattering rate has a 
crossover with respect to u as 

^ f e _7r/9 for 770 > u) 
\ log for rjQ<^u 

Perturbation theory with respect to disorder can be applied to the regime with 770 <C w, at least for 
quantities like the conductivity. The correlation length is not affected qualitatively by this cross over, as 
one can see from the matrix elements of r 2 

($ u/2 |^|$_ u/2 )=-4a;- 2 ((%Vff 2 )/ fc E i (ifit-^i) for Vo » " , 

[l for 770 <C lu 

which provides an effective correlation length lo^ 1 for (<S>E jT .j<!>*—E,rj). As we go away from the Dirac 
point, the correlation length decreases. On the other hand, the correlation function of the local density 
of states D rr / clearly distinguishes both regimes, since the corresponding Fourier transform is 

[q) ~ \ DiS(q) for t? < u 
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An interesting consequence of the "universal" uj 2 behavior of the matrix element of r\ is a constant 
microwave conductivity in both regimes: 



This behavior must be seen in contrast to the conventional Drude behavior, where the real part of 
conductivity decays like uj~ 2 . It is a consequence of the correlated scattering between the upper and the 
lower part of the Dirac cones (i.e., it is a manifestation of the zitterbewegung) . This effect should be 
experimentally observable, since it survives also in the presence of disorder. 

Our results can be summarized by the statement that properties of graphene at the Dirac point are 
ruled by long-range correlations. At the distance lo from the Dirac point they decay on the scale uj^ 1 . 
This has physical consequences, like a constant microwave conductivity, which also has been observed 
experimentally 13J . Disorder of variance g creates a characteristic scale 770 = e~ 7r / 9 , which separates 
the behavior in a vicinity lu <C 770 of the Dirac point from another one away from the Dirac point with 
u> ^s> 770 • Standard perturbation theories and simple approximations can be applied to the latter because 
then the behavior is ruled by one-particle properties. 
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Appendix A: Conductivity 

The conductivity per site on a lattice with N sites can be evaluated within the Kubo formalism and gives 
(cf. Eq. (3) in Ref. [25]) 

Re(a kk ) = -^-n [ Tr [[H, r k ]5(H -E + u)[H, r k ]5(H - E)] f ^ E + ^ ~ f ^ E) dE 
IV fx J UJ 

= _£L^ 2 J Tr M(H _ E + u/2)r k S(H -E- u /2)] M E + »M ~M E ~ »M dE 

with the Fermi function at inverse temperature (3: fp(E) = 1/(1 + exp(0E)). Now we consider 

Tr [r k 8(H - E + uj/2)r k S(H - E - u/2)] = ^ r k r' k Tr 2 [8 rr ,(H - E + uj/2)8 r , r {H - E - u/2)] 

We can write 

2r k r l k = -{r k -r' k f + rl+r' k 2 
Moreover, S(H -E + u/2)5(H - E - u/2) = for u ^ 0. Thus we obtain 

r kr' k Tr 2 [6„>(H — E + u/2)5 r , r {H -E- u/2)] 
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\ $>k - r fe) 2Tr 2 [Srr>{H -E + w/2)5 r , r {H -E- uj/2)] 



r,r' 



This gives 



-TTU) 



I X> - r '^ Tr2 ^r< H - E + »/*)St>AH -E- uj/2)} ff>& + U>/2) -fp{E - U>/2) ^ 

f p (E) ~ e(-s) , 



27V« 

Since for low temperatures 
we get 



e 2 ^ r" /2 

Y\{r k -r' k ) 2 Tr 2 [5 rr ,(H - E + uj/2)5 r , r {H - E-uj/2))dE 

~^ J-u/2 



The Dirac delta functions can be expressed by one-particle Green's functions: 

S(H + uj/2) = -L[G(w/2 - ie) - G(uj/2 + ie)} 

ZIYI 

such that 

-Re(crfefc) ~ ^(w) + <7 H (uj) 



with 



■2 r-u/t 

^(uj) = —uj / V(r fe - r' k ) 2 Tr 2 [G rr >(uj/2 -E- ie)G r , r (-uj/2 - E - ie) 

4ft J-u/2 ^ 

+G rr ,(uj/2 - E + ie)G T > r {-u/2 - E + ie)]dE 



and 



2 ru/2 

a u (uj) = -—uj / V(r fc - r' k ) 2 Tr 2 [G rr , (u/2 -E- ie)G r , r {-uj/2 - E + ie) 

4/l J-u/2 Z 



e 2 r 12 

u u (uj) = -—uj / Y(r k - r k ) 2 Tr 2 [G rr ,(uj/2 -E- ie)G r , r (-u/2 - £ + ie) 

4/1 J-ul2„ 



a j2 

+G rr -(w/2 - £ + ie)G r > r {-u/2 - E - ie)]dE 
and because of G(z) = -<t 3 G(-z)(T3 is 

+G rr /(-u;/2 + E - ie)G rlr {uj/2 + E + ie)]dE . 



Appendix B: matrix elements of the energy states 

Starting from a (localized) state ^(0), we can allow the state to evolve in time: 

= e- lHt V(0) 

The question is how the states ^±e can be reached by this evolution, and how this is influenced by 
scattering due to disorder. The contribution of the state with energy ±E to the time evolution is 
obtained by the Fourier transformation t ±E of ^f(t) for positive time t > (since the wave function 
did not exist for t < 0): 

/>oo />oo 

$ ±E = / e (±lE ~ e)t ^(t)dt = / e (±tE - e)t e- tHt dW(0) 
Jo Jo 

= -i(H T E - iey^iO) = -iG(TE - ie)*(0) . 
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Since H is Hermitean (i.e. iP = H), the complex conjugate of the wave function is 

&* ±E = iG(TE - ie)***(0) = iG(^E + ie) T **(0) , 
where G T is the transposed of G. The matrix element of r 2 between $±e is 

lpE\^-E)=Y,^**ri*- E ri^ (53) 

r r y 

if we assume that ^E'(O) is localized at the origin of the lattice r — 0. In the presence of disorder, this 
expression should be averaged with respect to the latter: 

(* E \'l\*-E)=Y f 1 i(*E,r3*-E,ri) = E T l ^ Ji' (~ E ~ ( E + ie )> (°) f ■ ( 54 ) 

r,j r,j,j' 

It describes a correlation between the states in the upper and in the lower band, if they evolve from the 
same initial state ^(O). These results can be summarized to the relation 

E r '< Tr2 [Gr^-E - ie)G r(-E' + ie)]) = {* E \fi\*E>) 

r 

if |*/(0)| 2 = 1 for f = 1,2. 

Appendix C: Inverse participation ratios 

The inverse participation ratio is related to the fourth moment of the normalized eigenfunction ^ with 
eigenvalue Ep. as 

P (2) =<El*fc,r| 4 <K£-£fc)> • (55) 
k 

This expression vanishes in the dclocalizcd regime. The latter is plausible when we estimate the delocal- 
ized wavefunction by l^fc.rl 2 ~ 1/N on a lattice with N sites, taking into account normalization 

Ei^i 2 = i. 

r 

This implies 

E E \^kA 2 5{E -E k )~Y< S ( E - Ek ) ~ N P ' 

r k k 

where p is the spatially averaged density of states. Moreover, for P' 2 ) we get from Eq. (|55|) a vanishing 
expression for N ~ oo: 

E Ms(e - Ek )~l,Y, S ( E - E *) ~ £ • 

fe fe 
P^ 2 ) can also be related to D rr by using the relation [3U] 

Urn e V |1V| 2 <5 £ (P - P,.) V |*fc/, r | 2 « e (£? - E„) = V l*fc,r| 4 <5(P - P fe ) 

e ^0 * — ' * — ' £~ 

which implies 

P (2) = lim eD rr . 
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Appendix D: Random- walk expansion 



The presence of the Grassmann field in the action S' of Eq. ((29)) enables us to write the functional integral 
J exp(— S')T>[ip, ip] in terms of a dense system of self-avoiding random walks (or polymers), whose density 
is controlled by the symmetry-breaking term. Starting from S' we either use UU' of Eq. (|2"0|) or its 
approximation given in Eq. (|30[) . In both cases the structure of the action with symmetry-breaking 
parameter e is 

S' = Sq + e VvVv + J]] a rr 'ip r 'ip r i + b rr >ip r , ip r ifj r i'ip r > . 

So is a term without Grassmann field. The quartic term (which does not appear in the approximated 
field of Eq. (|30p) can be expressed by a quadratic term that couples to a random field. Thus the general 
structure in terms of the Grassmann field reads 

S' = S + e /J VvVv + ^ B rr >i/} r -iJ) r i , 

t r,r' 

where the coefficients B rr / are random and connect only nearest-neighbor sites (this is a consequence of 
the symmetry). The integration with respect to tp can be performed first and gives just a determinant 



J exp(-S')£>[V>] = e- So det(e + B) 



Then the determinant can be expanded with respect to B, which creates random walks, whereas e creates 
an environment of identical points. In other words, a lattice site r is either occupied by an e or it is visited 
by a random walk with jump rate B rr i between r and a nearest-neighbor site r' . The random walks do 
not intersect themselves or each other and must be closed. This is a consequence of the integral over the 
Grassmann field or equivalently it is determined by the structure of the determinant. The contributions 
of the random walks depend on the weight of its elements B rr i, relative to the value of e. For e ~ (i.e. 
near the Dirac point) the lattice is completely covered by closed self-avoiding random walks, which can 
also be considered as the regime of dense polymers. Next, we need to integrate with respect to tp. In 
general, this will affect the weight of the random walks. In the special case of S' in Eq. ([31) we obtain 

J exp(-S')V[ip,(p] = J e- Sa det{K{e - iE) + zK{0)z)V[ip] , (56) 

where z is a diagonal matrix with matrix elements z rr = (<pn. r + ip22,r)/v and 

So = K rr>(t - iE)(-Lpii^ r tpiiy + Lp22,r^22,r>) ■ 

r,r' 

The expansion of the determinant produces terms with products z l r r r = YlrK^ii.r + f22.r)/i r l)\ lr (lr — 
0, 1, 2), where l r > is created by the second term in its argument z rr K rr i(0)z r i r i . Integration with the 
weight factor e _So , where at least one site r appears with l r > 0, gives a vanishing Gaussian integral 



J e- So Y[(ip lhr + ^r) lr =0 



since the two quadratic terms in So appear with opposite signs (i.e. ipn must be integrated along the 
imaginary axis). Therefore, the second term in the determinant of Eq. (|56p does not contribute to the 
functional integral. 
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Figure 1: Scattering on the Dirac cone: schematic scattering process between energy states &e-u>/2 
and <&E+u/2 a s well as between $_ w /2 and < i ) w /2- This type of scattering is relevant for the microwave 
conductivity a 11 ^) in Eq. ([5]). 
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Figure 2: Scattering rate (full curve) and density of states at the Dirac point (dashed curve) for Dirac 
fermions. Both quantities practically vanish over a wide range of disorder g. 




Figure 3: Fourier transform of the local density of states correlation function, which diverges like q 
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